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Simulated results from a 2+lD relativistic viscous hydrodynamic model have been compared 
to the experimental data on the centrality dependence of invariant yield, elliptic flow (^2), and 
hexadecapole flow (114) as a function of transverse momentum (pr) of charged hadrons in Au-Au 
collisions at ^snn ~ 200 GeV. Results from two types of initial transverse energy density profile, 
one based on the Glauber model and other based on Color-Glass-Condensate (CGC) are presented. 
We observe no difference in the simulated results on the invariant yield of charged hadrons for the 
calculations with different initial conditions. The comparison to the experimental data on invariant 
yield of charged hadrons supports a shear viscosity to entropy density ratio (rj/s) between to 0.12 
for the 0-10% to 40-50% collision centralities. The simulated V2{pt) is found to be higher for a fluid 
with CGC based initial condition compared to Glauber based initial condition for a given collision 
centrality. Consequently the Glauber based calculations when compared to the experimental data 
requires a lower value of r/fs relative to CGC based calculations. In addition, a centrality dependence 
of the estimated rj/s is observed from the V2{pt) study. The V4{pt) for the collision centralities 0-10% 
to 40-50% supports a rj/s value between - 0.08 for a CGC based initial condition. While simulated 
results using the Glauber based initial condition for the ideal fluid evolution under estimates the 
V4{pt) for collision centralities 0-10% to 30-40%. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 

Heavy-ion collisions at the Relativistic Heavy Ion 
Collider (RHIC) have provided evidence for the for- 
mation of a hot and dense QCD matter [l-iSi]. This 
presents an unique opportunity to study the trans- 
port properties, like shear viscosity to entropy den- 
sity ratio (ij/s), of the QCD matter. There are two 
main theoretical approaches to estimate the value of 
rj/s from the experimental data. One based on a mi- 
croscopic approach as in transport theory [6l-[l0| and 
other related to a macroscopic approach through rel- 
ativistic viscous hydrodynamic calculations jlll - [l8| . 
In this vifork, we will compare our results from a rel- 
ativistic 2-fl dimension viscous hydrodynamics to 
recent high statistics experimental data on elliptic 
(^2) and hexadecapole (U4) flow of charged hadrons 
measured by the PHENIX CoUaboration [H 
The experimental observables related to azimuthal 
anisotropic flow are found to be sensitive to shear 
viscous effects. The shear viscosity decreases the 
anisotropy of the fluid velocity. Hence V2 and W4 
as a function of transverse momentum (px) are ex- 
pected to decrease with the increase in the value of 
i]/s. 

One of the main uncertainties in the estimation of 
r]/s using a viscous hydrodynamics simulation is due 



to the choice of the initial conditions [13, [13) [HI ■ In 
this work we have considered two models, Glauber 
and Color Glass Condensate (CGC), to obtain the 
initial transverse energy density profile. For this 
study we have considered a smooth initial condi- 
tion, which does not vary event-by-event. Previ- 
ous work have shown that both the spatial and mo- 
mentum anisotropy are expected to be larger for a 
CGC based initial condition compared to Glauber 
model [17]. Hence for other similar conditions in the 
simulations, the calculations with CGC based initial 
condition is expected to give higher values of V2 com- 
pared to initialization based on a Glauber model. 
Earlier comparisons of viscous hydrodynamic sim- 
ulations with both CGC and Glauber initial condi- 
tions to the experimental data at RHIC can be found 
in Refs [13 [ll, [H [13 • In Ref [13, [3 , the experi- 
mental data used for comparison are the centrality 
dependence of multiplicity, (pr), Pt integrated U2, 
and minimum bias V2 vs. pr for charged hadrons 
in Au-Au coUisions at ^/snn = 200 GeV. In general 
it was observed that calculations with CGC based 
initial condition prefers a higher value oi rj/s com- 
pared to calculations with a Glauber based initial 
condition. In Ref [2lj, [l^j the authors have tried to 
explain the centrality dependence of V2 divided by 
the eccentricity with a viscous hydrodynamic model 
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for the QGP phase coupled to a transport model 
for the hadronic phase. Comparison of the experi- 
mental data to the calculations done for CGC initial 
condition supports a 77/s value ^ 0.16 - 0.24. While 
the corresponding comparisons for a Glauber model 
based initial condition supports a lower value of 77/s 
- 0.08 - 0.16. 

In the current work we compare the results from 
the viscous hydrodynamics simulations with two dif- 
ferent initial conditions (Glauber and CGC) to re- 
cent high statistics measurements of V2{pt) and 
Viipr) of charged hadrons in Au-Au collisions at 
y^SNN = 200 GeV for a broad range in collision cen- 
trality (from 0-10% to 40-50%) We also com- 

pare the simulated results to the measured charged 
particle invariant yields as a function of px for vari- 
ous collision centralities [20j . 

The paper is organized as follows. In the next 
section we discuss the formalism of viscous hydro- 
dynamic model used in this work. This includes a 
brief discussion on the energy-momentum conserva- 
tion and relaxation equations for shear stress. We 
present a detailed discussion on the initial conditions 
used in the calculations. The equation of state used 
and the freeze-out conditions are also presented. In 
section III we present a comparative study between 
calculations with Glauber and CGC initial condi- 
tions of various observables in the simulation. These 
includes the temporal evolution of shear stress, av- 
erage transverse velocity, and eccentricity. Section 
IV presents the comparison of viscous hydrodynamic 
simulations with different input values of 77/s for 
both Glauber and CGC based initial conditions to 
the experimental data on invariant yield versus px, 
V2{pt), and V4{pt) for various collision centralities. 
Finally in section V we present a summary of the 
work. 



II. VISCOUS HYDRODYNAMIC 
SIMULATION 

In a relativistic viscous hydrodynamics scenario, 
there are two-fold corrections to the ideal fluid hy- 
drodynamics. In presence of the dissipative pro- 
cesses, the energy momentum tensor contains ad- 
ditional dissipative corrections. The equilibrium 
freezeout distribution function used in the Cooper- 
Frey freezeout prescription [2^ also gets modi- 
fied. The first order dissipative correction to the 
energy- momentum tensor leads to acausal behav- 
ior 24]. The second order causal viscous hydro- 
dynamics due to Israel-Stewart is one of the most 
commonly used theory [2^. For the simulation re- 
sults presented here, we will follow the Israel-Stewart 
formalism for the evolution of a viscous fluid using 
the 2+ ID viscous hydrodynamic code "AZHYDRO- 



KOLKATA" [ii,!!^. Shear viscosity is the only dis- 
sipative process considered in our present study. We 
assume a net-baryon free plasma is formed in Au-Au 
collisions at midrapidity at ^/snn — 200 GeV. 

A. Conservation and relaxation equations 

The energy- momentum conservation equation and 
relaxation equation for shear viscosity in Israel- 
Stewart formalism is expressed as 

= 0, (1) 

- Ktt''^ + u''tt>'^]Dux. (2) 

Equation [1] is the conservation equation for the 
energy-momentum tensor, T^^ = (e -|- p)u'^u'^ — 
Pg^'^ + T^^'^ ■ £, P, and u are the energy density, 
pressure, and fluid velocity respectively, tt'^'^ is the 
shear stress tensor. Equation [5] is the relaxation 
equation for the tt^" . D — u^^d^ is the convec- 
tive time derivative, V<'"u''> = 1(^1*" + Vii^) - 
\{d^u^^){g^^'' — u^^u") is a symmetric traceless ten- 
sor. ?7 is the shear viscosity and is the cor- 
responding relaxation time. Assuming longitudi- 
nal boost-invariance, the above equations are solved 
with the code 'AZHYDRO-KOLKATA' in (r = 
— z^, X, y, -qs — l^'^f^f) coordinates. Where r 
is the longitudinal proper time, {t, x, y, z) are space- 
time coordinates, and rjs is the space time rapidity. 

B. Initial conditions 

The initial conditions used here includes the ini- 
tial energy density profile in the transverse plane 
(e(a;, 7/)), the initial time (t^), the transverse veloc- 
ity profile {vx{x,y),Vy{x,y)), shear stresses in the 
transverse plane {Tr^'^{x,y)) at r^. The value is 
taken as 0.6 fm. The 77/s values are also inputs to 
the viscous hydrodynamics simulations. We have 
taken the following temperature independent values 
for this work, 77/s = 0, 0.08, 0.12, 0.16, and 0.18. 

We have considered two different models for the 
calculation of initial energy density profile in the 
transverse plane. One is based on a two compo- 
nent Glauber model. At an impact parameter b, 
the transverse energy density is obtained from the 
following two component form 

e(b,2;, y) = eo[{l-Xh) (b, x, y)+XhNcoii (b, x, y)] 

(3) 
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where NpartO^, y) and Ncoiii^i, x, y) are the trans- 
verse profile of participant numbers and binary col- 
lision numbers respectively, ep corresponds to the 
central energy density in 6 = and does not de- 
pend on the impact parameter of the collision. The 
parameter Xh is the hard scattering fraction. Both 
eo and Xh are fixed to reproduce the experimen- 
tal charged hadron multiplicity density at midra- 
pidity. The iVpar-t(b, x, y) and Ncou{h,x,y) values 
are obtained using an optical Glauber model calcu- 
lation [isj. The value of xu is found to be 0.9 and 
the values of cq for various input values of rj/ s are 
given in the Table ID 



The other model commonly used to obtain ini- 
tial conditions for hydrodynamics is the Color-Glass- 
Condensate (CGC) approach, based on the ideas 
of gluon saturation at high energies [1^ [sO]. We 
have used the KLN (Kharzeev-Levin-Nardi) kx- 
factorization approach [3l|, due to Drescher et al. 

We follow references [13, and consider that 
the initial energy density can be obtained from the 
gluon number density through the thermodynamic 
relation, 



e(T,,XT,6) = C X 



-,4/3 



(XT,&) 



(4) 



TABLE I. Values of eo used in Glauber model and nor- 
malization constant C used in CGC model for initial 
transverse energy density. 



ri/s 


eo (GeV/Zm^), Glauber 


C (GeV/Zmi/^) qqq 


0.0 


43.4 


0.11 


0.08 


36.5 


0.095 


0.12 


32.5 


0.085 


Ifi 


27,7 


0.070 


0.18 


25.4 


0.065 





where ^a^^^y is the gluon number density evaluated 
at central rapidity Y = Q and the overall normaliza- 
tion C is a free parameter. C is fixed to reproduce 
the experimentally measured charged particle multi- 
plicity density at midrapidity. The values of C used 
in the simulations for different input values of rj/s 
are given in Table HI The number density of gluons 
produced in a collision of two nuclei with mass num- 
ber A is given by 



-HkT)V4;XT) (t>A(x2,{-pT 



kT)V4;XT), (5) 



where pt and Y are the transverse momentum 
and rapidity of the produced gluons, respectively. 
xi^2 = Pt ^ exp(±y)/-y/s is the momentum fraction 
of the colliding gluon ladders with ^/s the center of 
mass collision energy and ag^kx) is the strong cou- 
pling constant at momentum scale k^ = [k-xl- A/" is 
the normalization constant. The unintegrated gluon 
distribution functions are taken as 



(a;,fcj,;xT) 



1 



as{Q^) max((32^ fcf,) 



,^P{^t){1~x)\ 



(6) 

P(xt) is the probability of finding at least one 
nucleon at transverse position x^ and is defined 



ctTa ' 



where Ta is the 



as P{'^t) = 1 — (l J wnere i 

thickness function and a is the nucleon-nucleon 
cross section taken as 42 mb. The saturation 
scale at a given momentum fraction x and trans- 
verse coordinate xt is given by Q'^{x,x.t) 

icr)/P(xr) ' 
.53/fm^ 

is taken to be A = 0.28. 

Shear stresses tt^'^ is initialized to their corre- 
sponding Navier-Stokes estimates for the boost in- 
variance velocity profile, tt^^ = tt*'*' = 2ri/3Ti, tt^^ = 
[s^l- We have used r^r = Srj/ip (where rj and p 



2GeY' (^4^0-)) (Mi)\ The growth speed 



are the shear viscous coefficient and pressure) in our 
simulation, which corresponds to the kinetic theory 
estimates of relaxation time for shear viscous stress 
for a relativistic Boltzmann gas ^25j . The initial val- 
ues of Vx{x, y) and Vy{x, y) are taken to be zero. 

C. Equation of state 

In the present simulations we have used an equa- 
tion of state with cross-over transition at temper- 
ature Tc = 175 MeV The low temperature 
phase of the EoS is modeled by hadronic resonance 
gas, containing all the resonances with M^es <2.5 
GeV. The high temperature phase is a parametriza- 
tion of the recent lattice QCD calculation [SQ]- En- 
tropy density of the two phases are joined at T = 
Tc = 175 MeV by a smooth step like function. The 
thermodynamic variables pressure (p), energy den- 
sity (e), entropy density (s) etc. are then calculated 
by using the standard thermodynamic relations 



p(T)= r s{T')dT' 
Jo 

e{T)^Ts {T)-p{T). 



(7) 

(8) 
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D. Freeze-out condition 

The hydrodynamic expansion of the hot and dense 
matter leads to coohng of the system. After some 
time the mean-free path of the constituent becomes 
large/comparable to the system size. The system 
can no longer maintain the local thermal equilibrium 
and the momentum distribution of the particles re- 
mains unchanged after that. This is called freezeout. 
We use the Cooper-Frey algorithm at the freezeout 
to calculate invariant yields of the hadrons 23] . The 
freezeout temperature which is a free parameter in 
the hydrodynamics simulation is taken as T/=130 
MeV. The effect of different choices of freeze-out 
temperature on charged hadron px spectra and el- 
liptic flow is discussed in appendix \K\ 

As we have already pointed out there are twofold 
correction to the ideal fluid in the presence of vis- 
cous effects. The freezeout distribution function for 
a system slightly away from local thermal equilib- 
rium can be approximated as [1^ 

fneq {x,p) = feq (x, p) [1 -f- p)] , (9) 
I 



where 4>{x,p) << 1 is the corresponding deviation 
from the equilibrium distribution function feq(x,p). 
The non-equilibrium correction (pix^p) can be ap- 
proximated in Grad's 14 moment method by a 
quadratic function of the four momentum p^ in the 
following way [13, Issj 

0(x,p) =£-£^p^ + e,,,pV, (10) 

where £, e^, and e^i, are functions of p'^, metric ten- 
sor (7'"^, and thermodynamic variables. 

For our study where only shear stresses are con- 
sidered, <j){x,p) has the following form 

0(a:,p) =e^,p>^ (11) 

where 

^'"^ " 2(e + p)r2 ^^^^ 

As expected, the correction factor increases with in- 
creasing values of shear stress tThv at freezeout. The 
correction term also depends on the particle momen- 
tum. The Cooper-Frey formula [2^ for a non equi- 
librium system is 



where g is the degeneracy of the particle considered 
and c?Ep is the normal to the elemental freeze-out 
hypersurface. 

III. GLAUBER VERSUS CGC INITIAL 
CONDITION 

A. Space-time evolution 

Figure [T] shows the constant temperature contours 
corresponding to = 175 MeV and T/ = 130 MeV 
in the t-x plane (at ?; = 0) indicating the bound- 
aries for the QGP and hadronic phases respectively. 
The results are from the viscous hydrodynamic sim- 
ulations for Au-Au collisions at impact parameter 
7.4 fm and rj/s = 0.08. The solid red curves cor- 
responds to initial transverse energy density profile 
based on CGC model and the dashed black curve 
corresponds to results based on Glauber model ini- 
tial conditions. We observe that the lifetime of QGP 
and hadronic phases are slightly larger for the sim- 
ulations based on CGC initial conditions compared 
to Glauber based initial conditions. The spatial ex- 
tent of the hadronic phase is slightly smaller for the 



simulations with CGC initial conditions relative to 
Glauber based conditions. 



B. Temporal evolution of shear stress 

In presence of shear viscosity the thermodynamic 
pressure is modified. The tracelessness of shear 
stress tensor tt^"^ , along with the assumption of lon- 
gitudinal boost invariance ensures that at the initial 
time TT^^ and tt^^ components of shear viscous stress 
are positive. Consequently in viscous fluid the ef- 
fective pressure in the transverse direction is larger 
compared to the ideal fluid, for the same thermody- 
namic condition. It is then important to have some 
idea how various components of shear viscous stress 
tt'^^ evolves in space-time. We have considered tt^^, 
TT^y, and TT^y as the three independent components 
of shear stress n^'^ . This choice is not unique. 

The temporal evolution of spatially averaged tt^^, 
n^y, and tt^^ are shown in Fig.[2]for CGC (sohd red 
curve) and Glauber (black dashed curve) initializa- 
tion of energy density. All the three components of 
tt'^" becomes zero after time '--^ 7 fm irrespective of 
the CGC or Glauber model initialization. At initial 
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FIG. 1. (Color online) Constant temperature contours 
denoting the space time boundaries of the QGP and 
hadronic phases from a 2+ ID viscous hydrodynamic sim- 
ulation with ri/s = 0.08 for Au-Au collisions at impact 
parameter 7.4 fm. The quark-hadron transition temper- 
ature in the simulation is around 175 MeV and freeze- 
out temperature is taken as 130 MeV. The solid red 
curves are simulations with initial transverse energy den- 
sity profile based on CGC model while the dashed black 
curves correspond to initial conditions based on Glauber 
model. 
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FIG. 2. (Color online) The spatially averaged shear vis- 
cous stresses tt^^, tt''^, and tt^^ as a function of evolution 
time for Au-Au collisions at impact parameter 7.4 fm and 
r]/s — 0.08. The solid red and black dashed curves cor- 
responds to simulations with CGC and Glauber based 
initial conditions respectively. 



time the values of spatially averaged tt^^ and tt^^ 
are observed to be larger for CGC compared to the 
Glauber initialization. However, the difference van- 
ishes quickly ~ 3 fm. For tt^^ a noticeable difference 
is seen for CGC and Glauber model initialization 
within time ^ 6 fm. 
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FIG. 3. (Color online) Time evolution of spatially aver- 
aged transverse velocity {{vt))- The results are from a 
2+lD viscous hydrodynamic simulation with i]/s — 0.08. 
The solid red curve corresponds to simulated result with 
CGC based initial transverse energy density profile. The 
black dashed line is the simulated result with Glauber 
based initial conditions. 



C. Average transverse velocity and 
Eccentricity 

Figure [3] shows the temporal evolution of the spa- 
tially averaged transverse velocity (((wt))) of the 
fluid with Glauber based and CGC based initial 
transverse energy density profile with viscous hydro- 
dynamic simulations carried out for rj/s = 0.08. The 
simulation is done for Au-Au collisions at impact pa- 
rameter, b — 7.4 fm. The space averaged transverse 

velocity is defined as ((wt)) = ((7)) " ' '^'^^^^ 
^ , The angular bracket ((...)) implies 
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an energy density weighted average. Solid red curve 
is for CGC based initial condition and the dashed 
black curve is for the Glauber based initial condi- 
tion. We observe almost no change in the ((wt)) as 
a function of time for the two initial conditions stud- 
ied. This effect should be reflected in the slope of 
the invariant yield of the charged hadrons as a func- 
tion of transverse momentum being same for both 
the initial conditions. These results are discussed in 
section IV A. 

Figure H] shows the temporal evolution of the 
spatial eccentricity {Sx) and the momentum space 
anisotropy (Sp) of the viscous fluid (77/ s = 0.08) with 
Glauber and CGC based initial conditions for Au- 
Au collisions at impact parameter, b = 7.4 fm. The 
Sx which is a measure of the spatial deformation of 
the fireball from spherical shape is defined as 



(13) 
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FIG. 4. (Color online) The temporal evolution of spatial 
(ex) and momentum (gp) eccentricity for Au-Au colli- 
sions at b=7.4 fm. The solid red curves corresponds to 
viscous hydrodynamics (77/s = 0.08) simulated results 
with CGC based initial condition and the black dashed 
lines corresponds to results with Glauber based initial 
condition. 



and the £p which is a measure of the asymmetry of 
fireball in momentum space is defined as 



~ Jdxdy{T^ 



Tvv) ' 



(14) 



where T^^ and T^^ are the components of energy- 
momentum tensor T^". Solid red curve is for CGC 
based initial condition and the dashed black curve 
is for the Glauber based initial condition. We find 
both Ex and Sp are higher for the simulated results 
with CGC based initial condition compared to initial 
condition based on Glauber model. As the simulated 
elliptic flow V2 in hydrodynamic model is directly 
related to the temporal evolution of the momentum 
anisotropy, we expect the V2 for the CGC based ini- 
tial condition to be larger than the corresponding 
values for the Glauber based initial condition. These 
results are discussed in section IV B. 



IV. COMPARISON TO EXPERIMENTAL 
DATA 

The experimental data used for comparison to our 
simulated results are from the PHENIX collabora- 
tion at RHIC [H, [131 ■ The observables used are in- 
variant yield of charged hadrons, elliptic flow, and 
hexadecapole flow as a function of px for Au-Au 
collisions at pseudorapidity | ?y |< 0.35 for -^snn 
= 200 GeV. The high statistics recent PHENIX 
measurement of elliptic (fc — 1) and hexadecapole 
{k — 2) flow 19] are obtained using the formula 
V2k — {cos{2k{4>—^2))) after correction of the event 
plane resolution. Where (p is the azimuthal angle 
of the charged hadrons and '^2 is the second order 
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FIG. 5. (Color online) Invariant yield of charged hadrons 
as a function of transverse momentum at midrapidity 
for Au-Au collisions at ^snn = 200 GeV. The open cir- 
cles corresponds to experimental data measured by the 
PHENIX collaboration [23|. The lines represent results 
from a 2-1- ID relativistic viscous hydrodynamic model 
with a Glauber based initial transverse energy density 
profile and different rj/s values. The results are shown 
for five different coUision centralities 0-10%, 10-20%, 20- 
30%, 30-40%, and 40-50%. 



event plane constructed using event plane detectors 
in 1.0 <| 77 |< 3.9. The rapidity gap between the 
detectors used to measure the V2k and ^'2 ensures 
absence of significant A77 dependent non-flow corre- 
lations, which are also absent in our hydrodynamic 
simulations. ^2 for our simulation is along x axis. 
We compare below our simulated results on invariant 
yield, elliptic, and hexadecapole flow for five differ- 
ent collision centralities with input rj/s varying be- 
tween 0.0 to 0.18 to the corresponding experimental 
data. 



A. Invariant yield 

Figure [S] shows invariant yield of charged hadrons 
as a function of transverse momentum at midrapid- 
ity for Au-Au collisions at ^snn = 200 GeV for 
five different collision centralities (0-10%, 10-20%, 
20-30%, 30-40%, and 40-50%). The open circles 
are the experimental data from the PHENIX col- 
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FIG. 6. (Color online) Same as Fig.[5]but the 2+lD rela- 
tivistic viscous hydrodynamic simulations are done with 
a CGC based initial transverse energy density profile. 



laboration .20J. The simulated results are from the 
2+lD relativistic viscous hydrodynamic model with 
a Glauber based initial transverse energy density 
profile. The black solid, orange long dashed, pur- 
ple dash-dotted, magenta short dashed and green 
dotted lines corresponds to calculations with rj/s = 
0.0, 0.08, 0.12, 0.16, and 0.18 respectively. We find 
the 0-10% experimental data is best explained by 
simulation with rj/s — 0.0. Whereas data for colli- 
sion centralities between 20-30% to 40-50% supports 
a rj/s value within 0.08 to 0.12. 

Figure E] shows the same results as in Fig. [5] but 
the simulated results corresponds to 2-l-lD viscous 
hydrodynamic calculations with a CGC based ini- 
tial transverse energy density profile. The conclu- 
sions regarding the comparison between simulated 
results and experimental data are similar to that ob- 
tained for Fig. [5] This also means that the invariant 
yield of charged hadrons are not very sensitive to 
the choice of a Glauber based or CGC based initial 
conditions. The average transverse velocity at the 
freeze-out which determines the slope of the pT spec- 
tra was observed to be similar for the fluid evolution 
with Glauber and CGC based initial conditions (see 
Fig.©. 



FIG. 7. (Color online) Elliptic flow of charged hadrons 
as a function of transverse momentum at midrapidity 
for Au-Au collisions at ^/sm^ ~ 200 GeV. The open cir- 
cles corresponds to experimental data measured by the 
PHENIX collaboration [l^. The lines represent results 
from a 2-1- ID relativistic viscous hydrodynamic model 
with a Glauber based initial transverse energy density 
profile and different rj/s values. 



B. Elliptic flow 

Figure [7] shows the elliptic flow (^2) as a function 
of transverse momentum (pt) for charged hadrons 
at midrapidity in Au-Au collisions at y/swN = 200 
GeV. The results are shown for five different colli- 
sion centralities (0-10%, 10-20%, 20-30%, 30-40%, 
and 40-50%). The open circles are the experimental 
data from the PHENIX collaboration The sim- 
ulated results are from the 2+lD relativistic viscous 
hydrodynamic model with a Glauber based initial 
transverse energy density profile. The black solid, 
orange long dashed, purple dash-dotted, magenta 
short dashed, and green dotted lines corresponds to 
calculations with tj/s = 0.0, 0.08, 0.12, 0.16, and 0.18 
respectively. We find the experimental data prefers 
higher values of rj/s as we go from central to pe- 
ripheral collisions. While 0-10% collision centrality 
experimental data is best described by ideal fluid 
(ry/s = 0.0) simulation results, those corresponding 
to 40-50% collision centrality is closest to simulated 
results with rj/s ~ 0.18. 



8 




FIG. 8. (Color online) Same as Fig.[7]but the 2+lD rela- 
tivistic viscous hydrodynamic simulations are done with 
a CGC based initial transverse energy density profile. 
Also shown for comparison is the result with Glauber 
based initial conditions for ideal fluid evolution. 



Figure |S] shows the same results as in Fig. [7] but 
the simulated results corresponds to 2+ ID viscous 
hydrodynamic calculations with a CGC based ini- 
tial transverse energy density profile. Also shown for 
comparison the simulated results for ideal fluid evo- 
lution with Glauber based initial conditions. We find 
the V2{pt) for CGC based initial condition is larger 
compared to corresponding results from Glauber 
based initial conditions. This can be understood 
from the fact that CGC based initial condition leads 
to a higher value of momentum anisotropy compared 
to Glauber based initial condition (as seen in Fig.S]). 
The general conclusion that the experimental data 
prefers a higher value of rj/s as we go from central 
to peripheral collisions as seen for viscous hydrody- 
namic simulations with Glauber based initial condi- 
tions also holds for those with the CGC based ini- 
tial conditions. However, we find from the compar- 
ison of experimental data to simulations based on 
CGC initial conditions that the V2{pt) data for 0- 
10% collisions is best explained for simulated results 
with rj/s between 0.08-0.12. This is in contrast to 
what we saw from the comparisons of data to simu- 
lations with Glauber based initial conditions, where 
the data preferred rj/s — 0.0 (see Fig. [7]). For more 



FIG. 9. (Color online) Hexadecapole flow of charged 
hadrons as a function of transverse momentum at midra- 
pidity for Au-Au collisions at y^SNN ~ 200 GeV. The 
open circles corresponds to ex peri mental data measured 
by the PHENIX collaboration . The curves represent 
results from a 2-(-lD relativistic viscous hydrodynamic 
model with both Glauber based and CGC based initial 
transverse energy density profile and different rj/s values. 



peripheral collisions (centralities beyond 20-30%), it 
seems data would prefer a higher value of 7//s ^ 0.18. 
We do not present simulation results for rj/s > 0.18 
as the viscous hydrodynamic simulated spectra dis- 
tributions show a large deviations from ideal fluid 
simulation results (see appendix |B|) . This leads to 
a breakdown of the simulation frame work which is 
designed to be valid for case of small deviations of 
observables from ideal fluid simulations. 



C. Hexadecapole flow 

Figure |9] shows the hexadecapole flow {V4) as a 
function of transverse momentum (pt) for charged 
hadrons at midrapidity in Au-Au collisions at ^snn 
— 200 GeV. The results are shown for five differ- 
ent collision centralities (0-10%, 10-20%, 20-30%, 30- 
40%, and 40-50%). The open circles are the exper- 
imental data from the PHENIX collaboration [l9|. 
Simulated results for only ideal fluid evolution us- 
ing Glauber based initial condition are shown (solid 
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black curve). While for the CGC based initial con- 
ditions the simulated results are shown for 77/s = 0.0 
(purple solid thick curve) and 0.08 (orange dashed 
curve). We do not present simulated V4 results for 
other 77/s values as these are much lower compared 
to the data. We find that W4(pt) from ideal hydrody- 
naniic simulations with Glauber based initial condi- 
tions under predict the experimental data for all col- 
lision centralities studied except for the most periph- 
eral collisions (40-50%) presented. This is in sharp 
contrast to the observation for V2{pt) (see Fig. [T]) 
under similar conditions. Comparsion between sim- 
ulated results with CGC based initial condition and 
experimental data shows that the preferred 77/s lies 
between 0.0 and 0.08 for the collision centralities 
studied. The rj/s values supported by the data on 
V2 and W4 using the simulated results presented here 
appears to be different. In this study we have used 
smooth initial conditions, a more realistic approach 
is to use a fluctuating initial condition and carry out 
event-by-event hydrodynamics. This will enable us 
to study the odd flow harmonics ^3 along with the 
even harmonics V2 and W4. The simultaneous de- 
scription of all these experimentally measured flow 
harmonics in a viscous hydrodynamics framework 
will probably provide a better estimation of 77/s. 

V. SUMMARY 

We have carried out a 2-t-lD relativistic viscous 
hydrodynamic simulation with two different initial 
conditions (Glauber and CGC) for the transverse 
energy density profile in Au-Au collisions at y^SNN 
= 200 GeV. The simulations are carried out for rj/s 
values between 0.0 to 0.18, using a lattice -|- hadron 
resonance gas model based equation of state which 
has a cross over temperature for the quark-hadron 
transition at 175 MeV. The shear viscous corrections 
are considered both in the evolution equations and 
freeze-out distribution function. 

We find that the temporal dependence of the aver- 
age transverse velocity of the viscous fluid is similar 
for both the initial conditions studied. The compo- 
nents of shear viscous stress are observed to have 
higher values for the simulations with CGC initial 
conditions compared to those for Glauber model ini- 
tialization at early times of fluid evolution (< 6 fm). 
The simulated invariant yield of charged particles 
as a function of transverse momentum is also found 
to be similar for the Glauber and CGC based initial 
conditions. The spatial eccentricity and the momen- 
tum anisotropy have larger values for simulations 
with CGC based initial condition compared to the 
corresponding values for Glauber based initial con- 
dition. The simulated elliptic flow is observed to be 
higher for calculations with CGC based initial con- 



ditions relative to those with Glauber based initial 
conditions, for a given collision centrality. 

We have compared our simulated results to the 
experimental data at midrapidity on the centrality 
dependence of invariant yield, V2, and W4 as a func- 
tion of pt of charged hadrons measured in Au-Au 
collisions at y^snn — 200 GeV. From the compar- 
ison to the Pt spectra of charged particles we ob- 
serve that the data supports a 77/s value between to 
0.12 for the 0-10% to 40-50% collision centralities for 
both the initial conditions considered. The V2{pt) 
experimental data requires a lower value of 77/s for 
simulations with Glauber model initialization com- 
pared to the CGC based initial conditions. For both 
the models of initial conditions the V2{pt) data in- 
dicates a centrality dependence in the estimated 77/s 
value, with peripheral collisions preferring larger val- 
ues. The experimental data on V4{pt) for the col- 
lision centralities 0-10% to 40-50% supports a 7;/s 
value between - 0.08 for a CGC based initial con- 
dition. While simulated results using the Glauber 
based initial condition for the ideal fluid evolution 
under estimates the V4{pt) for collision centralities 
0-10% to 30-40%. Simulations with Glauber model 
initial conditions explain the V4{pt) data for 40-50% 
collisions with 77/s = 0.0. The observation associated 
with V4 is different from V2, with V4 data preferring 
smaller values of 77/s. 

There are further scopes of improvement on the 
simulated results presented here. Recent experimen- 
tal measurements of odd and higher order azimuthal 
anisotropic flow [39l - l4l| suggests that a fluctuating 
initial condition needs to be considered. It is ex- 
pected that the input 77/s to the hydrodynamic sim- 
ulations has a temperature dep endence in both the 
QGP and hadronic phases [ll|. Although large un- 
certainties still exist in the QCD computations of 
7;/s for the QGP phase. A more precise estimation 
of 7;/s would require the viscous fluid simulations to 
also consider bulk viscosity and vorticity effects [i^] . 
Both of which are expected to be non-zero for the 
system formed in high energy heavy-ion collisions 
and may affect the observables like V2- A proper 
prescription for bulk viscous freeze-out correction is 
still under debate in literature [13, |H, , while im- 
plementation of vorticity in viscous hydrodynamic 
simulations have just started to be investigated. We 
plan to consider some of these effects in the near 
future. 
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FIG. 10. (Color online) (a) Invariant yield and (b) el- 
liptic flow of charged hadrons as a function of transverse 
momentum at midrapidity for Au-Au collisions at ^/sm^ 
— 200 GeV. The open circles corresponds to experimen- 
tal data measured by the PHENIX collaboration [l9l ]. 
The lines represent results from a 24- ID relativistic vis- 
cous hydrodynamic model with a Glauber based initial 
transverse energy density profile with rj/s — 0.08 and 
different T/ values. 



Appendix A: Freezeout temperature 

We have studied the effect of different freeze-out 
temperature on the charged hadron invariant yield 
and elliptic flow. Figure [TU] shows the invariant yield 
(panel - a) and elliptic flow (panel - b) of charged 
hadrons for 20-30% centrality Au-Au collisions at 
^snn = 200 GeV. All the simulated results are for 
rj/s — 0.08 using the same Glauber based initial con- 
dition but with three different freeze-out tempera- 
tures, Tf — 110 (red dashed curve), 130 (black solid 
curve), and 150 MeV (blue dash-dotted curve). 

The slope of the pr spectra increases as Tf de- 
creases. This is because of higher radial velocity 
gained due to longer duration of evolution of the 
system for the lower freeze-out temperature. The ex- 
perimentally measured pt spectra is best explained 
for simulations with input parameters as specified in 
section II and Tf = 130 MeV. 

For the current simulations the effect of different 
Tf values studied is observed to be small on ellip- 
tic flow of charged hadrons. This could possibly 
be due to saturation of the value of the momentum 
anisotropy at the early time of evolution. 



Appendix B: Viscous correction 

There are two kinds of dissipative correction to 
the ideal fluid simulation. First the energy momen- 



tum tensor contains a viscous correction and the 
freezeout distribution function is also modified in 
presence of the dissipative processes. The viscous 
hydrodynamics model is applicable when the dis- 
sipative correction (both in the energy-momentum 
tensor and freezeout distribution function) is small 
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FIG. 11. Dissipative correction to the invariant yield of 
the charged hadron as a function of pt for 0-10% (black 
solid curve) and 40-50% centrality (black dashed curve) 
Au-Au collisions at -^snn = 200 GeV. The results are 
shown for r;/s=0.18 and Glauber based initial condition. 



compared to the corresponding equilibrium value. It 
is then implied that the relative viscous correction 
{5N/Neq) is small for the Israel-Stewart's hydrody- 
namics to be applicable, where Neq is the invari- 
ant yield for system in local thermal equilibrium, 
and 6JSI — cpprp \ viscous (pp-p \ equiiibrium • In figure 
[11] the relative viscous correction SN/N^q for ry/s = 
0.18 are shown for 0-10% (black solid curve) and 40- 
50% (black dashed curve) collision centrality. We 
observe that the relative shear viscous correction is 
quite large (~ 50%) at pt ^l-G GeV for 40-50% cen- 
trality (dashed curve in fIgure[TT|) . The solid curve in 
the same figure shows the relative viscous correction 
to the invariant yield of charged hadron for Au-Au 
collision for 0-10% centrality at ^snn = 200 GeV. 
A higher value oi rj/s will introduce a larger viscous 
correction and eventually the viscous hydrodynam- 
ics framework will no longer be applicable. 
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